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ABSTRACT 

We report on our analysis of HST/NICMOS snapshot high resolution im- 
ages of 255 stars in 201 systems within ~10 parsecs of the Sun. Photometry 
was obtained through filters F110W, F180M, F207M, and F222M using NIC- 
MOS Camera 2. These filters were selected to permit clear identification of 
cool brown dwarfs through methane contrast imaging. With a plate scale of 
76 mas/pixel, NICMOS can easily resolve binaries with sub-arcsecond separa- 
tions in the 19 / .'5xl9 / .'5 field of view. We previously reported five companions to 
nearby M and L dwarfs from this search. No new companions were discovered 
during the second phase of data analysis presented here, confirming that stel- 
lar /substellar binaries are rare. We establish magnitude and separation limits 
for which companions can be ruled out for each star in the sample, and then per- 
form a comprehensive sensitivity and completeness analysis for the subsample of 
138 M dwarfs in 126 systems. We calculate a multiplicity fraction of 0.01qq% for 
L companions to M dwarfs in the separation range of 5 to 70 AU, and 2.3~^qj% 
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for L and T companions to M dwarfs in the separation range of 10 to 70 AU. We 
also discuss trends in the color-magnitude diagrams using various color combina- 
tions and present astrometry for 19 multiple systems in our sample. Considering 
these results and results from several other studies, we argue that the so-called 
"brown dwarf desert" extends to binary systems with low mass primaries and is 
largely independent of primary mass, mass ratio, and separations. While focus- 
ing on companion properties, we discuss how the qualitative agreement between 
observed companion mass functions and initial mass functions suggests that the 
paucity of brown dwarfs in either population may be due to a common cause and 
not due to binary formation mechanisms. 

Subject headings: binaries: close — infrared: stars — stars: brown dwarfs, low 
mass, solar neighborhood, statistics 



1. Introduction 



The mass function, multiplicity fraction, and the mass-luminosity relation are three of 
the most important characteristics of a stellar or substellar population. However, all three 
remain poorly constrained for Very Low Mass (VLM) stars. Although the lowest mass stars, 
the M dw arfs, dominate the Galaxy in numbers and comprise the majority of our stellar 
neighbors (IHenry et al.ll2006[ ). not a single M dwarf is visible to the naked eye. Over the last 
two decades, advances in observational astronomy ha ye made a thorough study of these faint 



stars possible. Empirical mass-luminosity relations ( IHenry fc McCarthy! Il993l ; IHenry et al. 



19991 : iDelfosse et al.ll2000f ) have achieved a high degree of reliability for early to mid M dwarfs 
with progress continuing for later M dwarfs at the end of the main sequence. Large sky sur- 
veys such as the 2 Micron All Sky Survey (2MASS) in the ne ar infrared (ISkrutskie et al. 



2006 ) and the Sloan Digital Sky Survey (SDSS) in the optical ( York et al. 200ol ) have pro 



vided a wealth of new data for population studies, but lack the angular resolution necessary 
to investigate the Multiplicity Fraction (MF) and Companion Mass Function (CMF) at 
separations corresponding to the short periods necessary for determining dynamical masses. 



The discovery of GJ 229B, the first unequivocal brown dwarf (jNakajima et al.lll995l ). 
followed by hundreds of othersj, raised fundamental questions about our understanding of 
low mass star formation. Are VLM stars and brown dwarfs products of a single mecha- 
nism of (sub) stellar formation applicable to a wide range of masses? Or, do brown dwarfs 



A current list of known L and T dwarfs is maintained at www.dwarfarchives.org. 
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constitute a fundamentally different population? Does the trend in the stellar mass func- 
tion producing more stars at lower masses continue into the brown dwarf regime, there- 
fore making them more numerous than their stellar cousins? What do multiplicity proper- 
ties, such as the overall multiplicity fraction and the separation distribution, tell us about 
the environments in which VLM stars and brown dwarfs were born? These are some of 
the fundamental questio ns that have only rece ntly been addressed through a combina- 
tion of sky su rveys (e.g., 



searches 



Allen et al. 



.g. jReid fc Gi^islll997l: 



2008 



Bochanski et al. 
20071: 



2010) , wi de separation common proper motion 



A 



Close et al 



en 



20031 : 



fc Reidl l2008h . high resolution multiplicity surveys 



Gizis et al. 2003 



Lowrance et al. 2005; Reid et al. 



Metchev &: Hillenbrand 20091) and the establishment of trig onometric parallaxes for a 



large sample of objects (e.g., iDahn et al.ll2002l ; iHenry et al.l 120061 ) 



Thorough characterization of any stellar population requires the study of a volume 
limited sample. In an effort to better understand these properties, we have conducted an 
HST/NICMOS snapshot program imaging 255 objects in 201 star systems with trigonometric 



parallaxes placing them within ~10 pc of the Sun 



methane imaging ( iRosenthal et al.lll996t iTinney et al. 



Tabl e 1). We used the technique of 



20051 ) to clearly distinguish cool brown 



dwarf companions. In 2004 we repo rted the detection of four M dwarf companions and one 
binary L dwarf in a triple system (jGolimowski et al.l l2004bl . hereafter G04). With small 
infrared contrasts ranging from ~0 magnitudes (GJ 1001BC) to 4.5 magnitudes (GJ 84AB), 
the companions we reported in 2004 were relatively bright. We have since carried out a deeper 
search of the data, establishing formal sensitivity limits for the detection of companions in 
the field of each primary target and extending the limiting magnitude differences routinely to 
11 at separations of 3'.'0, 8 at l'.'O, 4 at 0'.'4, and 2 at 0'.'2 (§5.3, Figure 6a). Having completed 
the deeper search of the data with no further detections, we now report on the magnitude 
and separation limits to which we can rule out companions for each object in our sample. We 
also discuss what the lack of additional brown dwarf detections tells us about the multiplicity 
fraction of systems with VLM secondaries in these mass and separation regimes. 

We describe the general characteristics of our sample and discuss how the observed 
sample relates to our current knowledge of the solar neighborhood in §2. Instrumental 
aspects of the observations relevant to obtaining our sensitivity limits are briefly reviewed 
in §3, and a detailed discussion of our Point Spread Function (PSF) insertion method for 
testing the sensitivity of the search is given in §4. We discuss photometric trends in our 
color system and note several benchmark objects in §5.1. We report new astrometric data 
for 19 known binary systems in §5.2, and discuss the sensitivity of the search in §5.3. In §5.4 
we establish a sub-sample of 126 M dwarf systems for which we calculate the multiplicity 
fraction, including substellar companions, from 5 to 70 AU based on companion detections 
(or lack thereof) and completeness arguments. We discuss what our results mean in the 
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context of the "brown dwarf desert" in §6.1 and §6.2. Finally, we end by comparing young 
cluster multiplicity studies, estimates of the Galactic disk luminosity function, and our results 
in §6.3, and summarize our conclusions in §7. 



2. General Characteristics of the Sample 

Our target list was designed to provide a sample that is representative of the solar 
neighborhood. Because this survey was an HST snapshot program, our targets were pulled 
from a large pool of selected targets in order to fill smallgaps in HST's observing schedule. 
The REsearch Consortium on Nearby Stars (RECONSJo is engaged in an effort to obtain 
a census and thorough characterization of the population o f stars within 25 p c of the Sun, 



with a particular concentration on stars closer than 10 pc (IHenry et al.l 120061 ) . In order to 
be a member of the RECONS 10 pc sample, an object must have a trigonometric parallax 
greater than 100 mas, with an error smaller than 10 mas. We used the RECONS 10 pc 
sample as a starting list for our search and allowed the HST snapshot scheduling process to 
effectively select a random subsample from the 10 pc sample. Table 1 summarizes several 
tallies of the observed sample. These observations comprise 69% of the RECONS 10 pc 
sample (epoch 2012.0), including main sequence stars, white dwarfs, L and T dwarfs, but 
excluding extrasolar planets. We note that because trigonometric parallaxes for nearby stars 
are constantly being updated, 17 objects in 13 systems originally included in our search are 
no longer members of the 10 pc sample. We still include their data as individual stars in this 
paper, but exclude them from statistical considerations in order to keep the sample volume 
limited. 

Figure 1 shows the spectral type distribution of the NICMOS snapshot sample. Out of 
the 218 resolved objects within 10 pc we observed, 138, or 63%, are M dwarfs. This number 
is a very close match to the M dwarf fraction in the RECONS 10 pc sample, which is 248 
out of 357 objects, or 69% (epoch 2012.0). The preponderance of M dwarfs in our sample 
means that even though the sample is a random representation of the solar neighborhood, 
it focuses on the spectral type that is least scrutinized by RV companion searches and open 
cluster imaging searches. By studying nearby M dwarfs, which comprise a disk rather than 
cluster population, we are mapping the brown dwarf desert in a largely unexplored region. 



information about RECONS and periodically updated versions of the 10 pc census are available at 
www.recons.org. 
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3. Observations and Data Reduction 

G04 describe technical aspects of the observations in detail. We give a brief summary 
here and highlight the aspects that are most relevant in achieving the sensitivities we later 
quote for each individual target. 

We obtained direct images of each target using NICMOS Camera 2 (NIC2) through 
four near infrared filters during cycles 7 (1997-1998) and 11 (2002-2003). NIC2 has a 
plate scale of 0"076 pixel" 1 and a field of view of 19 r /5xl9'/5 (Viana et al. 2009, Thatte et 
al. 2009)|. Targets were imaged through the F110W, F180M, F207M, and F222M filters, 
centered at 1.10/iin, 1.80/iin, 2.07/nn, and 2.22/mi, respectively. Because HS T observations 
are not subject to atmospheric absorption, the NICMOS filters are not defined to sample 
atmospheric transmission windows in the way that ground based near infrared filters are. 
The resulting filter set is non-standard when compared to ground based systems, but allows 
the user to construct a color scheme that is more suitable for the underlying physics being 
investigated. Figure 2 shows the transmission curves for the selected NICMOS filters, with 
the 2MASS J, H, and K$ filters also plotted for comparison. The four filters in this survey 
were selected to detect the strong CH 4 absorption bands observed in T dwarf spectra at 
1.7/zm and 2.2^xm, in effect imaging in and out of these absorption bands. Depending on 
the filter choice, there is a drastic color shift of up to three magnitudes for T dwarfs. A 
late T dwarf appears blue in F110W—F180M (0.0 to —1.0, Figure 3a) whereas it is red in 
F180M—F207M (1.0 to 2.0, Figure 3c). Because no background source is likely to have such 
a strong color shift, T dwarfs are readily identified in this color scheme. This technique 
is commonly known as methane im aging; and has been used to successfully i dentify brown 



dwarfs in photometric observations ([Rosenthal et al.lll996t iTinney et al.ll2005l ) 



By centering the targets on the detector, we searched for companions within a radius of 
9'.'5, except for a small (~1" in diameter) artifact due to the coronagraphic hole on the upper 
left quadrant of the deteetofl A few targets had large coordinate uncertainties, in most cases 
due to poorly constrained high proper motions. These targets were not properly centered in 
the field of view, and are specified in the notes to Table 2. Although some of our primary 
targets are very bright (e.g. Sirius, Vega, Procyon), we did not use the coronagraph because 
it would make the acquisition process too long for a snapshot program and its peripheral 



3 HST/NICMOS documentation, including the NICMOS Instrument Handbook and the NICMOS Data 
Handbook, is currently available from the Space Telescope Science Institute at www.stsci.edu/hst/nicmos 

^HST's roll orientation during a given exposure is constrained by the need to keep the solar arrays facing 
the Sun. Consequently, the position angle of the coronagraphic hole with respect to celestial north, as well 
as the image's overall orientation, varies widely among the images of our targets. 
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position in the detector would severely limit our search radius. Placing the primary target 
behind the coronagraph would also add uncertainty to the measurement of the position angle 
and separation of any binary systems. Even with saturated central targets, we could still 
search for companions, albeit with a lower sensitivity closer to the central target (Table 2). 

We coadded two sets of exposures for each target, resulting in a total exposure time of 
64s for the Fl WW and F180M filters and 128s for the F207M and F222M filters. Saturation 
of bright targets and cosmic ray hits were minimized by using NICMOS's multi-accumulate 
(MULTIACCUM) mode, which reads the detector in a non-destructive manner at predeter- 
mined time intervals. In the event of saturation or a cosmic ray hit, the NICMOS pipeline 
scales the value from unaffected readouts so as to obtain the approximate value due to the 
astronomical source. Only pixels that saturate or are hit by a cosmic ray before the first 
readout at 0.303 s are lost. For targets that were bright enough to saturate during the first 
readout, we obtained photometry by using PSF fits. 

Because there is no background atmospheric glow, the extended PSF of the primary 
target is the dominant source of background flux obscuring any fainter objects in the field 
of view. We subtracted a properly scaled PSF of another star of similar spectral type 
and brightness from the survey fro m the PSF o f each target. A detailed discussion of the 



PSF subtraction process is given in iKrist et al.l (119981 hereafter K98). The quality of the 



subtraction varied from target to target and depends primarily on whether or not a good PSF 
match could be obtained. The PSF varies with target color, telescope focus, and the position 
of the NICMOS cold mask (K98). We were always able to find an isolated star whose PSF 
was used as the reference for PSF subtraction. If the PSF reference had been a close binary 
or if it had been contaminated by background sources, we would have noticed a physically 
unrealistic negative PSF in the subtracted image. We then performed aperture photometry 
on the primary target as well as any other sources in the field of view using standard IRAF 
routines and the aperture corrections for encircled energy fraction listed in Table 2 of K98. To 
verify the validity of the aperture corrections, we performed the photometry of the crowded 
field of LHS 288 (31 sources, Figure 4), varying the aperture from three to six pixels (0'.'23 
to C/46). The photometry agreed to <0.03 magnitudes in all bands, regardless of aperture. 
For the final photometry we chose a six pixel aperture, except in cases when a crowded field 
or a source near the edge of the field required a smaller aperture. 



4. Determining the Sensitivity of the Search 



We define the "sensitivity" of the search as the extent to which we can detect or rule 
out the existence of a companion to a given star at a given separation and image contrast 
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Am. The sensitivity varies from target to target and is influenced by the overall brightness 
of the primary target, the quality of the PSF subtraction, the image filter, intrinsic detector 
noise, and the prominence of detector artifacts. For each image these factors interact in a 
complex way, thus making it difficult to draw generalizations about instrumental sensitivity 
for the survey as a whole. We have therefore devised a method to measure the sensitivity 
achieved for each target at various separations, and quote individual results in Table 2. 

Because HST is not subject to atmospheric effects, its images are inher ently stable, 



thus facilitating PSF modeling. We used Tiny Tim 6.3 (IKrist fc Hook! 119971 ) to simulate 
generic NIC2 stellar PSFs through the four filters used in the search. The properly scaled 
model PSFs were inserted into the PSF subtracted images of the primary targets to test 
our ability to detect companions at a range of contrasts and separations using a customized 
IDL code. At sub-arcsecond separations, we inserted a single companion at separations of 
0'.'2, 0'.'4, 0'.'6, and 0'.'8 and a varying range of contrasts at random position angles (Figure 
5a). The PSF insertion code automatically excluded the strong diffraction spikes present 
in well-exposed NICMOS images at 45°, 135°, 225°, and 315°. The residual flux from the 
PSF subtracted primary target was set to zero at a radius interior to the position of the 
artificially inserted companion to facilitate visual inspection. At separations of l'/O or greater, 
we produced an image where artificial companions were arranged in a radial pattern around 
the PSF subtracted primary (Figure 5b). This pattern tested the sensitivity at separations 
of l'/O, 2'.'0, 3'.'0, and 4'.'0 at contrasts typically incremented from 6 to 13 magnitudes. The 
simulated images and their surface plots were then visually inspected. In both regimes, an 
artificially inserted companion was considered detectable if it was visible in the simulated 
image and if a surface plot around the companion indicated that the artificial PSF retained its 
characteristic stellar shape, with its peak clearly above the background noise, corresponding 
to a typical signal-to-noise of 3—5. Although automating the PSF recovery process (e.g., 
by using a cross-correlation algorithm) would have saved a considerable amount of time, we 
were not convinced that automated methods would appropriately distinguish between real 
astronomical sources and residuals of the central star's PSF subtraction, which can at times 
mimic star-like profiles. 



5. Results 

Other than the five companions reported in G04, which focused on individual discoveries, 
we detected no further new companions during this second phase of our search. We now 
report on the photometry, astrometry, and search sensitivities attained during the survey. 
With a statistically robust sample of 255 stars surveyed, these results allow us to make 
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assessments of the multiplicity of stars in the solar neighborhood from a stellar population 
perspective. 



5.1. Color-Magnitude Diagrams 

We constructed color-magnitude diagrams for the twenty-four different color-magnitude 
combinations from our observations through the four filters. Because our sample includes 
only four certain substellar objects (GJ 1001 B and C, GJ 229B, and 2MASSI J0559191- 
140448), we used synthetic photometry from the spectra of known L and T dwarfs to better 
determine the form of the substellar sequence in this color space. These values were ob - 



tained using flux-calibrated, near infrare d spectra (IGeballe et al.l 120021 ; iKnapp et al.l 120041 ) 



weighted mean trigonometric parallaxes ( iGolimowski et al.l l2004bl and references therein) 



and the NICMOS Exposure Time Calculator produced by STScI. Figure 3 shows four color- 
magnitude diagrams that are particularly well suited for mapping the stellar and substellar 
main sequence. Main sequence targets and the thirteen white dwarfs in the survey are la- 
beled with large dots. In these diagrams, we initially assume that any object in the field of 
view of a primary target is a companion and therefore shares the primary's trigonometric 
parallax. If the assumption is correct, the object will fall within the stellar or substellar se- 
quence. Background objects, labeled with small dots, appear as having unrealistically faint 
absolute magnitudes and tend to cluster at the bottom of the diagrams. 

The trends in the F110W-F180M and the F110W-F222W colors (Figures 3a and 3b) 
clearly indicate that the onset of CH4 absorption happens sharply around the L6 spectral 
type, where the colors turn blue. Although any single diagram may show an overlap between 
the substellar sequence and the brighter background objects, the degeneracy is broken when 
we consider that L and T dwarfs follow different trends from the background sources in 
different color combinations. The most dramatic example of these shifts appears in Figures 
3c and 3d, where methane imaging causes a large shift from blue to red for the T dwarfs 
while the background sources show little change. 



5.1.1. Benchmark Objects 



GJ 1245 ABC (labels 1, 2, and 4 in Figure 3) is an interesting system containing three 
low mass components. In particular, GJ 1245C (4) is one of the latest M dwarf s for which a 
dynamical mass is known. With a mass of O.O74±O.O13M ( IHenry et al.lll999l ). this object 
lies close to the theoretical hydrogen burning mass limit. 



-9 - 



G 239-25B (label 3) was discovered during the first phase of this searc h, and the im- 



plicat ions of the multiplicity of the G239-25 system are discussed in G04. iForveille et al. 



( 120041 ) assign it a spectral type of L0±1 based on near infrared spectra. This spectral clas- 
sification makes G 239-25B an important benchmark of the M/L transition at the bottom 
of the main sequence. Its proximity in color space to GJ 1245C re-enforces the importance 
of both objects as benchmarks. 

GJ 1001BC (labels 5, 6, and 7) was resolved as a double L4.5 dwarf, and is discussed 
in detail in G04. The components of the system are nearly equal in luminosity, and we plot 
them both individually (6 and 7) and combined (5) to illustrate how an equal flux binary 
appears in the sequence. When compared to the L/T sequence outlined by the synthetic 
photometry, both components of GJ 1001BC lie just before the strong shift towards the blue 
that happens as a result of the onset of CH 4 absorption. Their positions at this turning point 
are most easily seen in the F110W—F222M color (Figure 3b). We are currently working to 
refine the parallax of GJ 1001 ABC, and to obtain dynamical masses for the BC pair. 

Finally, 2MASSI J0559191-140448 (T4.5, label 8) and GJ 229 B (T6, label 9) are the 
only T dwarfs imaged in the survey, and serve as confirmations that the sequence outlined 
by the synthetic photometry agrees with real photometry. Whereas GJ 229B is a companion 
to the M0.5V dwarf GJ 229A, 2MASSI J0559191-140448 is an isolated brown dwarf. Its 
positions in panels a and b of Figure 3 illustrate how a mid T dwarf can easily be mistaken 
for a white dwarf when more color combinations are not used to break the degeneracy. 



5.1.2. Background Objects with Companion- Like Colors 

Figure 3 shows that there are several sources having colors that mimic the colors of 
substellar companions in one or more panels. The ambiguity is often accentuated when 
analyzing data sets with simpler color combinations that were not designed a priori to dis- 
criminate substellar objects (e.g. 2MASS JHK S ). Interstellar reddening considerations are 
particularly useful in identifying false companions. Because the distance horizon of our 
search is only ~10 pc, any bona fide companions should not have appreciable reddening 
in the near infrared. Conversely, distant main sequence or giant stars may have significant 
reddening in the F110W-F180M, F110W-F207M, and F110W-F222M colors, which may 
place background objects in the color space occupied by L and T dwarfs. The degeneracy 
is broken when considering the F207M—F222M and especially the F180M-F207M colors, 
where the narrow spectral coverage reduces the reddening of distant main sequence sources 
(Figures 3c and 3d). Table 3 lists cases where the distinction between a background object 
and a putative companion was particularly subtle based on colors alone. The white dwarfs 
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as a group mimic late L and early T dwarfs in F110W-F180M and F207M-F222M, but 
the degeneracy is broken in F180M-F207M. 

5.2. Astrometry of Known Binaries 

High resolution images of nearby binary systems present opportunities to map relatively 
short period orbits and therefore determine dynamical masses. While actual orbital mapping 
is beyond the scope of this work, we report the astrometry for select systems in Table 4. 
In order to be listed in Table 4, both components of the system must have been imaged 
simultaneously in the same NIC2 field of view, and the centroids must be determined to a 
precision better than ±1 pixel. The values we report are the weighted averages of separations 
and position angles measured from the PSF centroids in all filters for which saturation did 
not prevent reliable centroiding. In the simplest case of non-saturated and non-overlapping 
PSF cores, we adopt a centroiding error of ±0.1 pixel (G04). Twelve out of the 19 pairs 
listed in Table 4 meet these criteria. The other seven pairs are either very closely separated 
stars for which the PSF cores overlap significantly or have central pixel saturation. In either 
case, the centroiding was determined using PSF fits. With the exception of the M dwarfs, 
the majority of binaries in our survey had their PSF cores saturated beyond the point where 
we could compute meaningful astrometry. The precise value of the NICMOS plate scale 
varied during HST cycle 7 (1997—1998) due to cryogen expansion that distorted the dewar 
housing the detectors. To calibrate the plate scales for our observations, we used the values 
tabulated by the Space Telescope Science Institute based on routine monitoring of crowded 
star fields. For separations, the errors listed in Table 4 take into account the four centroiding 
uncertainties (x a , y a , Xb, Ub) added in quadrature. For position angles, the errors take into 
account the propagated centroiding errors. 

5.3. Results from the Sensitivity Search 

Table 2 lists the faintest detectable absolute magnitudes for putative companions at a 
range of angular separations from each target star in the survey. The distances and spectral 
types listed are based on the best trigonometric parallaxes and spectral type estimates avail- 
able in the literature or unpublished trigonometric parallaxes recently measured or improved 
by our group. It is important to note that each line in Table 2 shows the results of one PSF 
insertion simulation, and does not necessarily correspond to a single star. As described in 
the notes column, a single PSF insertion simulation may have been done around the two 
components of a resolved system if their separation was small or if their contrast was large 
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enough for the primary to dominate the field. Because of these situations, the number of 
entries in Table 2 is not meant to add up to the sample counts in Table 1. The reader is 
referred to Table 1 for overall statistics of the sample and to Table 2 for data on individual 
targets. 

All targets were inspected for real companions visually in all four bands over the entire 
field of view. Several factors must be considered when choosing the best filter for the PSF 
insertion simulations. Out of the four filters used in the search, the Fl WW and F180M filters 
are the most suitable for close separations (< Of! 4) due to their narrower PSFs when compared 
to the F207M and F222M filters. Whereas L dwarfs are brighter in F180M than in F110W, 
the T dwarfs are much fainter in F180M due to methane absorption. Although the F110W 
band produces the narrowest PSFs due to its shorter wavelength and is an intrinsically bright 
band for T dwarfs, we chose to report the sensitivities in the F180M band for two reasons. 
First, our uniform exposure time scheme (§3) causes brighter targets to saturate out to 
several pixels in the FllOWband even in 0.303s, decreasing our ability to probe the smallest 
separations. Second, the width of the F110W PSF is comparable to the NIC2 pixel scale, 
causing a sharp spike on the central pixel (K98, Table 2). In low signal-to-noise situations 
it becomes difficult to distinguish the F110W PSF from bad pixels or other sharp artifacts 
introduced during the PSF subtraction process. As discussed in §5.4.2 our sensitivity limit 
falls mostly in the L dwarf regime for sub-arcsecond separations, and in the T dwarf regime 
for wider separations. Based on comparisons in particularly clear images, we estimate that 
using F110W instead of F180M would increase our sensitivity by ~1 magnitude, but would 
pose an unacceptable risk of false detections at close separations. We therefore uniformly 
report sensitivities for all separations in F180M, but emphasize that those values can be 
safely transformed to F110W limiting magnitudes for separation greater than l'.'O by adding 
1.0 magnitude to the F180M limits in Table 2. Because late T dwarfs appear the faintest in 
F180M, a detection in that band also implies detection in the other three bands, therefore 
providing the color information needed to characterize the object. Listing our simulation 
results in the F180M band therefore maximizes the instrumental dynamic range of the images 
while still providing the sensitivity needed to characterize T dwarfs. 



5.4. The M Dwarfs 



Of the 188 star systems imaged within 10 pc in this survey, 126 systems have M dwarfs 
as the primary (or single) component^. Because these M dwarfs were selected randomly from 
a volume limited sample based on their trigonometric parallaxes, the sub-sample lends itself 
well to statistical considerations. We now apply the sensitivity limits in Table 2 to derive 
the multiplicity fraction for M dwarfs under several scenarios. 



5.4-1- Establishing Search Completeness for M Dwarfs 

Figures 6 a and b show the ranges in sensitivities obtained for M dwarfs at each of the 
eight angular separations probed by the PSF insertion simulations. In Figure 6b, we used 
the known distance to each target to convert contrasts into absolute magnitudes, and relate 
these absolute magnitudes to the spectral types of putative companions. Because sensitivity 
is a complex function of contrast, instrumental background, apparent magnitude, and the 
quality of the PSF subtraction, there is a significant spread about the mean values quoted 
in Figure 6. Overall, we would detect companions with AF180M=2.5 to 10.2 magnitudes at 
separations of 0'.'2 to 4'.'0, respectively. 

In order to transform our observational sensitivities (Figure 6a) to astrophysical param- 
eters, we substitute physical separations in AU in place of angular separations and apply the 
statistical relation between physical separation and semi-major axi s for a sample of binarie s 



with random inclinations and eccentricities, < a >= 1.26 < p > ( IFischer fc Marcyl Il992l ). 
obtaining Figure 7. The large plus signs in Figure 7 indicate the 90% detection limits for 
semi-major axes ranging from to 40 AU, binned in 2 AU increments. We assume a flat 
contrast curve for sensitivities beyond 40 AU. Because of the large factor in distance covered 
by this volume-limited search, the 90% detection limits in physical separation are effectively 
established by the most distant stars in the sample. It is possible to boost sensitivity at 
closer physical separations by establishing a closer distance horizon for the search, at the 
expense of overall sample size. We examined the effect of using a closer distance horizon 
for calculating sensitivity limits, and came to the conclusion that it is more important to 
maintain a robust sample, especially because more sensitive but much smaller studies have 



5 GJ 169. 1AB is an M4.0V/white dwarf binary. Although the brighter M4.0V component is generally 
considered to be the primary component, the current situation does not reflect the components' masses or 
spectral types at the time of stellar formation and main sequence evolution, when the current white dwarf 
was much more massive and luminous than the M dwarf. We therefore do not consider GJ 169. 1AB to be a 
system with an M dwarf primary. 
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already been done (e.g., IClose et al.l 120031 ) . We emphasize that Table 2 contains all data 
necessary for different statistical formulations, and is available in machine readable format 
in the online version of this paper. 

We also considered the effect that the small field of view of NIC2 (l$!5xl$!5) has 
on sample completeness at large physical separations. Figure 8 is a histogram displaying 
the number of all M dwarfs, including resolved system secondaries, sampled within 10 pc 
(N=141) as a function of outer search radius, binned in 10 AU increments. While all M 
dwarfs were probed to semi- major axes as close as 5 AU@, only the farthest 12 targets were 
probed at semi-major axes greater than 120 AU. In order to retain the statistical significance 
of our sample, we consider only physical separations corresponding to mean semi-major axes 
between 5 and 70 AU (100% to 79.4% complete), and divide the number of companions 
found in the bins from 40 to 70 AU by that bin's completeness fraction. 



5.4-2. The M Dwarf Multiplicity Fraction 

Table 5 lists companions to M dwarfs in our sample within our completeness range of 5 
to 70 AU that were re-detected in our search or are new companions discovered during this 
search and published in G04. Combining these known binaries to the null detections and 
sensitivity limits we present in Figures 6—8, we now present formal multiplicity fractions 
for three distinct combinations of companion types and ranges in semi-major axes. These 
results are summarized in Table 6. In each case, the la confidence intervals were calcu- 



lated using the binomial distribution approach outlined by iBurgasser et al.l (120031 ) . This 
approach is preferable to traditional Poisson statistics whenever the probability distribution 
is non-Gaussian. In each of our three different scenarios discussed below the multiplicity 
fraction is low enough that even with the sample of 126 systems, the probability distribution 
is not symmetric about the central peak value because proximity to the limiting case of a 
multiplicity fraction of zero causes a sharper drop-off towards the lower limit of the proba- 
bility distribution (Figure 9). Given a multiplicity fraction e m , the probability distribution 
of finding n binaries in a sample of iV systems is governed by 

AT' 

P{n) = UN • ., £(1 - e m ) N ~ n . 
n ] .(jy — n)\ 



6 As noted in Table 2, GJ 15A, LHS 224AB, GJ 623AB, and GJ 644ABD are M dwarfs for which core 
saturation prevented the establishment of a sensitivity limit at 0'.'2. A single M dwarf system, GJ 747AB, 
saturated out to 0'.'4. All of these cases correspond to statistically corrected semi-major axes smaller than 5 
AU. 
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This relationship can be inverted to solve for the probability distribution of a given multi- 
plicity fraction given the observational results N and n, yielding 

P'(e m ) = (N+l)P(n), 

which can then be integrated numerically to find the lower and upper limits of e m correspond- 
ing to 68% (la for a Gaussian distribution) of the area under the probability distribution 
curve, as shown by the shaded areas in Figure 9. 

The M Dwarf Multiplicity Fraction for MOV to M9V Companions at Separations of 5 to 70 

AU 

At an inner search radius of 5 AU, our search is 90% sensitive to Mfisom < 11.2, correspond- 
ing to early L spectral types (Figure 7a). Eleven of the 12 known companions listed in Table 
5 are M dwarfs meeting this sensitivity criterion. Five of these companions lie between 40 
and 70 AU, where the completeness of the search is reduced due to the limited field of view. 
Placing these five systems into the separation bins shown in Figure 8 yields 2 systems in the 
40-50 AU bin, 2 systems in the 50-60 AU bin, and one additional system in the 60-70 AU 
bin. Dividing these numbers by the fractional completeness of these bins (0.957, 0.879, and 
0.794) and summing the results yields 5.62. We then transform the multiplicity obtained at 
90% confidence level to a true volume limited multiplicity fraction by dividing 11.62 (the 
sum of 5.62 and the remaining companions from 5-40 AU) 

by 0.9, obtaining 12.91. Rounding this number up to 13, we see that we would likely have 
recovered 2 additional real companions with separations ranging from 5 to 70 AU. Applying 
the binomial distribution, we conclude that the multiplicity fraction for M dwarf companions 
orbiting M dwarf primaries at semi-major axes from 5 to 70 AU is e m = 10. 3^^% (Figure 
9a). 

The M Dwarf Multiplicity Fraction for L0 to L9 Companions at Separations of 5 to 70 AU 

Although our search did not detect any L dwarf companions within 10 pc and in the sep- 
aration regime of 5 to 70 AU0 it is possible to assign a multiplicity fraction based on 
completeness arguments. Figure 7a shows that at 5 AU, the detection rate for L dwarfs is 
only ~50%. It is not possible to obtain a truly volume limited multiplicity fraction in this 
separation range. We therefore constrain the sample to include only the 51 systems for which 



7 GJ 1001 B and C are beyond 10 pc (jHenrv et al 



200G). 
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the detection of an L9 companion at 5 AU is possible. Applying the binomial distribution, 
we obtain a multiplicity fraction of of e m = 0.0^^% (Figure 9b). An alternative approach 
is to maintain the volume limited nature of the sample by increasing the inner limit of the 
separation range. From Figure 7a, the inner radius at which >90% of the systems were 
probed is 12 AU. We therefore calculate a volume limited multiplicity fraction for L0 to L9 
companions to M dwarfs of e m = 0.0tlx,% valid at separations ranging from 12 to 70 AU. 



The M Dwarf Multiplicity Fraction for L0 to T9 Companions from 10 to 70 AU 

Our sensitivity to T dwarfs at close separations is diminished due to their intrinsic faintness. 
We therefore restrict the inner search radius to 10 AU, where the search was 90% sensitive 
to L dwarfs and ~50% sensitive to late T dwarfs. At separations beyond 12 AU, Figure 7b 
indicates considerable scatter in the 90% sensitivity limits. Based on the trend on Figure 7b, 
we adopt a 90% sensitivity limit of Mfuow= 17.5, corresponding to spectral type ~T9. One 
T6 dwarf, the class prototype GJ 229B, was detected at an inferred semi-major axis of 55.3 
AU. Following the same approach we used for the L dwarfs, we calculate the multiplicity 
fraction for a sub sample as well as for the volume limited sample. There were 43 systems for 
which a late T dwarf detection at 10 AU was possible. This sub-sample yields a multiplicity 
fraction of e m = 2. 3^;°% The complete sample is sensitive to late T dwarfs at separations ^ 
14 AU. We therefore calculate a volume limited multiplicity fraction of e m = 0.81q \% valid 
at separations ranging from 14 to 70 AU. 



6. Discussion 



6.1. Sensitivity to Companion Masses 



Estimating masses for field brown dwarfs is a difficult problem. W hereas the masses of 
main sequence stars can be estimated from mass-luminosity relations ( IHenry &: McCarthy 
19931 ; iHenry et al.lll999uDelfosse et al.ll2000l ). brown dwarfs are constantly cooling, and there- 
fore have a mass-luminosity-age relation. Such a relation has not yet been established empir- 
ically. Currently, the best way of estimating brown dwarf masses is by correlating spectral 
types to effective temperatures, and then checking the effective temperature against evolu- 
tionary model predictions, assuming a certain age for the brown dwarf in question. This 
approach is heavily model dependent, and the end result of such calculation can at best 
serve as a guideline for the mass range for a particular object. With this caveat in mind, we 
now apply this approach to the limiting spectral types we report in Table 6. 
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Assuming a mean age of 3 Gyr for the nearby L dwarf field population (ISeifahrt et al. 



20101 ) . the effective temperatures for brown dwarfs of spectral types L3, L 5, L8, T5, and T7 



are estimated to be 2000K 1750K, 1500K, 1200K, and 900K, resp ectively (IGolimowski et al. 



2004al ; ICushing et al.ll2008l ). Adopting the evolutionary models of IChabrier et al.l (120001 ). we 
estimate approximate masses of 0.073 M , 0.070 M , 0.057 M , 0.052 M , and 0.040 M 
for spectral types L3, L5, L8, T5, and T7 (Table 6). The last number has considerable 
uncertainty due to the steeper decline in effective temperatures for subtypes later than ~T5 
and the need to extrapolate the Chabrier models at low temperatures. We therefore adopt 
O.O4OM at 3 Gyr as a guideline for the minimum mass detectable by our search. We note 
that the scatter in age in the nearby field population is likely to cause a large dispersion in 
the masses of detectable objects. Unless there are further data indicative of the age of an 
individual brown dwarf, the mean value we adopt here should be used with extreme caution. 



6.2. A Current Map of the Brown Dwarf Desert 



The idea of the brown dwarf desert continues to evolve. The term was originally used 
to describe the fact that radial velocity surveys of solar analogs detect an abundance of 
extra-solar planets but rarely detect brown dwarfs, ev en though a brown dw arf's higher 
mass makes its detection easier. In their seminal work, iMarcy &: Butler! ( 120001 ) found that 
<1% of main sequence Sun-like stars harbor brown dwarfs. Several other studies have 
since then obtained similar results fo r different ranges in separation, primary mass, and 
system age. lOppenheimer et al.l (120011 ) conducted the first successful search for brown dwarf 
companions, discovering the T dwarf prototype GJ 229B. Their infrared coronagraphic search 
of stars within 8 pc detected a single substellar object, from which they c autiou sly imply a 
stellar /substellar multiplicity fraction of ~1%. McCarthy fc Zuckermanl (120041 ) used Keck 
coronagraphy to search 102 nearby field GKM stars at separations from 75 to 1200 AU. They 
found one brown dwarf companion, and report a binary fraction of l±l9fl We note that 
their result agrees well within statistical uncertainties to our results (Table 6), suggesting 
a wide de sert with no significa nt change in the substellar companion fraction from 10 to 
1200 AU. iLuhman et aD f|2005h used HST's Wide Field Planetary Camera 2 (WFPC2) to 
survey 150 members of the young cluster IC 348 (~2 Myr) at separations of 120—1600 AU. 
Of these stars, 85 were in the mass range 0. 08—0.5 Mm, approximately corresponding to 



the mass range for main sequence M dwarfs (IHenry fc McCarthy! 1 19931 ) . They found one 
possible substellar companion to a low mass star, but note that it is not possible to ascertain 



8 Using the binomial distribution treatment we adopt in this paper, 1 detection out of 102 observations is 
equivalent to a multiplicity fraction of 2 %. 
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companionship due to the wide separation of this system (~1400 AU). Based on this finding, 
Luhman et all derive an upper limit of 4% for the substellar companion fraction of low mass 
stars. This result is again in very good agreement with our results, suggesting that there is 
little evolution in the multiplicity fraction of low mass stars after the first few million years, 
and again suggesting no significant change in the substellar companion fraction beyond 10 
AU. Regarding single objects, Luhman et al. find that 14 out of 150 objects are likely 
substellar based on evolutionary models (IChabrier et al.ll2000l ). They note that the fact that 
they detect ten times more isolated stars than isolated brown dwarfs in IC 348 indicates 
that the brown dwarf desert may not be limi ted to the formation of compan ions, but may 
also extend to the formation of single objects. iMetchev fc Hillenbrand! (120091 ) used adaptive 
optics on Keck and Palomar to survey 266 Sun-like (F5— K5) stars, and infer a brown dwarf 
companion frequency of 3. for separations of 28 to 1590 AU. Finally, direct imaging 

searches for planetary companions would be capable of detecting brighter brown dwarfs. 
Masciadri et aD fj2005h used VLT/NACO to search 30 young (<200 Myr) GKM stars and 
found no brow n dwarf or planet ary companions at separations larger than 36 AU. In a 
similar fashion, iBiller et al.1 (120071 ) used VLT and MMT to search 45 young GKM field stars 
at separations of 20—40 AU, and also found no brown dwarfs. Due to smaller sample sizes, 
the last two studies do not add significant constraints to the brown dwarf desert, but their 
null detections are certainly in agreement with constraints set by the larger studies. 

The sum of these studies, along with the results we present in this paper, indicate a 
consistent image of a brown dwarf desert that is mostly invariant with respect to the mass 
of the primary star, and which is valid for a wide range of separations ranging from 5 AU 
to 1600 AU. Whether the search i s sensitive to substellar comp anions to Sun-like stars at 



intermediate to large separations ( IMetchev fc Hillenbrand! 120091 ) . substellar companions to 



or a mixture oi 


young stars with 


(Masciadri et al. 


2005 




Biller et al. 



20071 ) the detection rate is always consistent with a stellar-substellar binary fraction on the 
order of a few percent. 



6.3. Is the Desert Real? 



The multiplicity fraction of Sun-like stars is ~50 % ( jDuquennoy fc Mavoiil99ll : iRaghavan et al.l 

2010h. The multiplicity rate for stellar companions to M dwar fs at all separations is ~30% 



( IHenry fc McCarthylll990l ; lHenrylll99ll ; iFischer fc Marcylll992l ). Based on our results (Table 



2a limits. 
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6) and the companion searches we discuss in §6.2, it is clear that stellar companions out- 
number brown dwarf companions by a factor >10. Does this paucity of brown dwarfs 



however, constitute a "real desert"? A few studies (e.g., iMetchev &: Hillenbrand! 12009 



Grether &: Lineweaverl 120061 ) have suggested that the dear th of brown d warf companions 



is a natural consequence of a well behaved, Salpeter-like (jSalpeterl Il955l ) universal Com- 
panion Mass Function (CMF) that tends to lower multiplicities at lower mass ratios, and 
that a real brown dwarf desert would only exist if the observed number of brown dwarf 
companions is significa ntly lower than what a universal CMF would predict. In particular, 



Grether &: Lineweaverl note that the overlap of the planetary CMF and the stellar CMF 



reaches a minimum at ~0.03 M Q , causing the observed paucity of brown dwarf companions. 
In our search, we test the hypothesis of a universal CMF by focusing primarily on low mass 
stars. As shown in Table 5, the twelve M dwarf binaries we detected between 5 and 70 
AU have primary masses ranging from ~0.6 to ~0.1 M . Figure 10 is a plot of the masses 
of the primary and the secondary components of this sample. Figure 10 shows that the 
mass ratios of low mass binaries tend to increase (i.e. approach equal mass components) as 
masses approach the hydrogen burning limit, thus excluding the formation of brown dwarf 
secondaries. Our completeness analysis demonstrates that this trend is not an observational 
selection effect. Indeed, detecting companions with higher contrasts is easier for intrinsically 
fainter primary stars, so the selection effect works against the trend noted in Figure 10. 
Reconciling our observations with the idea of a universal CMF would require this function 
to be rather restricted in the sense that it would not be a function of mass ratio, or would 
only be valid for Sun-like stars. For any reasonably broad definition, we conclude that devi- 
ations from a universal CMF do exist in the brown dwarf regime. The brown dwarf desert is 
therefore a reality whether one defines it in terms of total numbers or in terms of a deviation 
from a trend. We advocate that the concept of a universal CMF is probably not a useful 
representation of Nature. 



6.4. More Evidence for A Discontinuity at the Hydrogen Burning Limit? 



VLM binaries have a strong tendency towards high (i.e. unity) mass ratios (e.g. lBurgasser et al 



20071 ). The effect has been demonstrated t o be an intr insic characteristic of VLM stars and 
brown dwarfs through Bayesian analysis (lAllenl 120071 ) . Our results (Figure 10) show that 
mass ratios tend to increase as stellar masses approach the hydrogen burning limit, with the 
strong onset of nearly equal mass duplicity happening somewhere between 0.2 and 0.1 M . 
Other studies have also suggested that the basic population properties of Initial Mass Func- 
tion (IMF), CMF, and the binary separation distribution all ap pear to change sign ificantly at 
a mass of ~0.1 M , slightly above the hydrogen burning limit. I Close et al.1 (120031 ) conducted 
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an adaptive optics search of 39 VLM objects with spectral types ranging from M8.0V to L0.5, 
and found a mass distribution similar to the one shown in Figure 10 (see their Table 3). They 
also probed smaller separations than our formal limit of 5 AU, and found that whereas higher 
mass stars have a separation distribution peaked at 30 AU ( IDuquennoy &: Mayorlll99ll ). VLM 
binaries have a sep aration distribution peaked at 4 AU. Also, Bayesian analysis of several 
studies (lAllenl 120071 ) demonstrates that VLM and b rown dwarf binaries with separations > 
20 AU are extremely rare. We note that lClose et al.l probed s ignificantly smaller separations 
than we did, but did not establish formal detection limits. iKraus et al.l (120051 ) conducted 
a search for VLM binaries in the Upper Scorpius OB association, and also found results 
consistent with a discontinuity in the separation distribution at a mass of 0.1M o . 



In an analysis of data from several open cluster studies, iThies &: Kroupal (120071 ) demon- 
strate that the observed mass distribution is incompatible with the existence of an IMF that 
is monotonic about the hydrogen burning limit. They note that because stellar formation 
and stellar ignition are in principle unrelated processes governed by different areas of physics, 
there is no reason to expect that the IMF discontinuity would be caused by the onset of hy- 
drogen burning. They therefore allow for an arbitrary overlap of the stellar and brown dwarf 
components of the IMF, thus allowing for a smooth turnover. In light of our companion 
mass distribution for low mass stars (Figure 10), new developments in the hydrodynamical 



simulations of star cluster formation ( Bate 



2009 



2011), and new observations of young stellar 



clusters (IKraus et al.l 120081 l2011t lEvans et al.ll2012l ). we re-examine the nature of the IMF 
discontinuity at masses close to the hydrogen burning limit. 

The details of the mass function for older field objects close to the hydrogen burning limit 
are difficult to quantify. The difficulty is mostly due to the lack of a robust volume limited 
census of L and T dwarfs based on trigonometric parallaxes or reliable distance estimates 
(errors <20%). For the M dwarfs, the situation is more clear. Our recent res ults from the 
REC ON S 10 pc cens u s indi cate a minimum M dwarf space density of 0.057pc -3 (IHenry et al. 



2006 



10 



Cruz et al.l (120071 ) find a space density of 4.9 x 10 



3 pc 3 for M dwarfs later than 



M7V and a lower limit of 3.8 x 10 pc for L dwarfs. Assuming that field age (~1— 5 Gyr) 
brown dwarfs with masses slightly below the hydrogen burning limit are predominately mid 
to late L dwarfs (§6.1, Table 7), and that stars o f spectral type M7V or later have masses 
<0.1 M ( IHenry et al.lll999l ; iDelfosse et al.l 120001 ). the ratio of objects with masses above 
and below 0.1 M is 6.9. The shape of the M dwarf distribution in the RECONS 10 pc 



10 0.059pc -3 for epoch 2012.0. See www.recons.org for the latest numbers and analysis. Comparison of the 
10 pc sample with the 5 pc sample indicates that the 10 pc M dwarf sample is ~70% complete. We note, 
however, that an analysis of the RECONS sensitivity limits indicates that the assumption of a representative 
M dwarf sample within 5 pc may be significantly biased by statistics of small numbers. 
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census corresponds broadly to the distribution of our NICMOS sample, (Figure 1), with the 
drop-off happening at around spectral type M6V, corresponding to ~0. 1 Mm- Ev en if the 



actual density for L dwarfs is a few times greater than the lower limit of I Cruz et al.l ( 120071 ). 
there is still a significant difference in the number of stars versus brown dwarfs. 

Could the onset of core hydrogen fusion cause the discontinuity in the IMF and the 
CMF via a radiative feedback mechanism? We caution that our understanding of stellar 
formation processes in this mass range is rather limited from a theoretical as well as an 
empirical basis, so any explanation is tentat ive at best. We speculat e that if the onset of 
core hydrogen burning at ages from 3—5 Myr ( IChabrier fc Baraffdll997l ) has a significant role 
in hindering accretion, the star formation process would produce a discontinuity at masses 
> 0.075 M Q . Although it has been generally accepted that protostars acquire the bulk of 
their mass during the first 1 Myr, observations sho w that a sizable fract i on of substellar 
objects continue to accrete for a much longer time. iJayawardhana et al.l ( 120031 ) find that 
40%— 60% of brown dwarfs in young star forming regions with ages up to ~10 Myr show 
infrared excesses consistent with accretion. There is also observational evidence that at 
least some high mass brown dwarfs und ergo phases of strong accretion compar able to the 
T Tauri phases of more massive stars (IBouy et al.l 120081 ; IComeron et al.l |2010| ). For the 
highest mass proto-brown dwarfs, late accretion may be enough to ignite hydrogen fusion, 
or to otherwise significantly change the manner in which the young object interacts with 
its environment. More observations and theoretical work are needed to confirm or discard 
this hypothesis, in particular with regards to testable predictions of accretion rates. Even 
if late accretion brings the total mass of a proto-brown dwarf above the hydrogen burning 
limit, we lack a clear understanding of how the onset of core hydrogen fusion would hinder 
accretion. At ages of a few Myr, the vast majority of an object's luminosity comes from 
the release of internal gravitational ene rgy, so the onset of hydro gen burning would have a 
negligible effect on overall luminosity (I Chabrier fc B ar affd 1 1 9 9 7( ) . We note, however, that 



hydrodynamical cluster collapse simulations (IBatd 120091 ) are in good agreement with the 
stellar IMF and stellar CMF, but overpr oduce the n umber of brown dwarfs unless radiative 
feedback is incorporated into the model (jBatdl201ll ). The last model produces a cluster of 
stars and brown dwarfs whose statistical properties are very similar to those of observed 
young clusters, suggesting that radiative feedback is indeed an important mechanism in 
brown dwarf formation. If the discontinuity in the CMF and the IMF at 0.1 M Q stands up to 
further observational scrutiny, a strong convergence of theoretical models and observational 
evidence will be needed to prove this or other hypotheses. 



7. Conclusions 



We conducted a large, volume-limited, high resolution search for substellar companions 
around nearby stars, with a particular emphasis on M dwarf systems. By evaluating the 
completeness and sensitivity of the search, we have established the multiplicity fractions for 
M dwarfs listed in Table 6. We find a multiplicity fraction of O.O+l'ffi, for L companions to M 
dwarfs with semi-major axis ranging from 5 to 70 AU. Including T dwarfs down to spectral 
type T9 and restricting the inner search radius to 10 AU yields a multiplicity fraction of 
2.3lj$%. These rates are far less than M dwarf pairs, for which we found a multiplicity 
fraction of 10.3l2'i% f° r separations of 5 to 70 AU. Based on these results, we summarize 
the substellar multiplicity fraction for M dwarfs as being on the order of a few percent or 
less. As discussed in §6.2, several other multiplicity studies have reached essentially the 
same conclusion regardless of primary mass, the separations probed, or the sample's age 
estimate. The emerging picture is that of a pervasive "brown dwarf desert", hinting to 
origins that are largely independent of a binary's formation mechanism. By specifically 
focusing on low mass primaries, our study has weakened the case for mass ratio dependence 
in the formation of substellar companions. We add ours to a list of several studies that 
indicate that the companion mass function is truncated at a mass ~0.1 M Q , slightly above 
the hydrogen burning mass limit (§6.4, Figure 10). While the primary focus of this work is 
characterizing the companion population, we also note in §6.4 that the results we obtain for 
the stellar /substellar multiplicity fraction are consistent with estimates for the population 
density of isolated brown dwarfs. The similarity suggests that mechanisms causing the 
observed paucity of brown dwarfs, both as companions and as isolated objects, may be 
intrinsic to the brown dwarf formation process. Recent results from hydrodynamic cluster 
collapse simulations as well as evidence for T Tauri like accretion at ages of a few Myr make 
radiative feedback from recently ignited very low mass stars a good candidate mechanism 
for truncating the IMF and the CMF at masses slightly above the hydrogen burning mass 
limit. 
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Fig. 1. — Spectral type distribution of the 239 targets within 10 pc in the search. The sample 
constitutes 69% of the RECONS 10 pc sample (epoch 2012.0). 63% of the targets are M 
dwarfs, which is in close agreement with the M dwarf distribution of the 10 pc sample, 69% 
(epoch 2012.0). 
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Fig. 2. — Transmission curves for the four NICMOS filters used in the survey. The 2MASS 
filters are plotted with dotted lines for comparison. Although no individual NICMOS filter is 
a close match to a ground based equivalent, together they cover nearly the same wavelength 
range of ground based near infrared color systems. 
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Fig. 3. — Selected color- magnitude diagrams designed to detect substellar companions. The 
large dots are the primary targets of the search, including 13 white dwarfs. The small dots 
are background objects. Synthetic photometry of L and T dwarfs, as well as one M9 dwarf, 
is plotted using a label for spectral type, with the precise dot position at the center of the 
label. In these diagrams, all objects within the field of view of a primary target are plotted 
assuming a common parallax (i.e. companionship). Only if the assumption is correct would 
the object fit in the stellar or substellar sequence. The benchmark objects discussed in §5.1.1 
are labeled as follows: (1) GJ 1245A, (2) GJ 1245B, (3) G 239-25B, (4) GJ 1245C, (5) GJ 
1001BC (combined), (6) GJ 1001B, (7) GJ 1001C, (8) 2MA 0559-1404, and (9) GJ 229B. 
Panels (a) and (b) illustrate the drastic color shift around spectral type L6 caused by the 
onset of CH4 absorption. The reduced effect of interstellar reddening on background objects 
displayed in panels (c) and (d), as well as the large shift from red to blue for substellar 
objects, make these bands particularly useful for methane imaging. 
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(b) 

Fig. 4.— Survey images for LHS 288 (M5.5V) using logarithmic scaling. The frames 
illustrate typical survey images both before (a) and after (b) PSF subtraction. The ghost- 
like coronagraphic artifact is visible in the upper left hand corner, particularly in the F222M 
images. The highly structured PSF of the primary target dominates the field before PSF 
subtraction. 



-29 - 




(b) 

Fig. 5. — Examples of sensitivity simulations around GJ 213 (M4.0V, F180M=6.68). (a) 
To test sub-arcsecond separations, a mosaic of PSF insertions with several separations and 
magnitudes is created. In this figure the rows represent separations of 0'.'2 (bottom), and 
(X'4 (top). The columns represent apparent F180M magnitudes of 9, 10, and 11 from left to 
right. The artificial companions are visible at all three magnitudes for Of 4 but only at the 
brightest magnitude for Of 2. (b) PSF insertions are laid out in a radial pattern to test the 
sensitivity at separations of l'.'O and greater. Apparent F180M magnitudes range from 12 to 
19 in increments of 1, with the rays for 18 and 19 not detectable in this case. Separations 
are l'.'O, 2'.'0, 3'.'0, and 4'.'0. In both simulations the residuals of the PSF subtraction are set 
to zero at a radius interior to the artificial companions to facilitate detection. A thorough 
inspection requires using surface and contour plots. 
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Fig. 6. — Search sensitivities for the eight angular separations tested by PSF insertion sim- 
ulations. In both panels, the two numbers next to each cluster of points are the mean and 
standard deviation for that separation, respectively, (a) The ability to detect a companion is 
primarily determined by the angular separation and the components' Am. This instrumental 
representation has a lower standard deviation, but does not probe fundamental astrophys- 
ical parameters, (b) Transforming Am into absolute magnitudes yields a range of possible 
companion types detectable at each angular separation. The absolute F180M magnitude for 
select spectral subtypes is taken from the synthetic photometry displayed in Figure 3b. 
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Fig. 7. — (a) Search sensitiv ity displayed as a func tion of absolute F180M magnitude and mean semi-major axis, 
assuming < a >= 1.26 < p > (IFischer fc Marcylll992l ). Each dot represents the sensitivity derived from a PSF insertion 
around an M dwarf (Table 2). The range of limiting absolute magnitudes is significantly wider at close separations 
because all targets were probed at close physical separations, whereas only targets close to our distance limit of 10 pc 
could be probed at wide physical separations given NIC2 's small field of view. Contrast is also more strongly dependent 
on overall brightness at close angular separations. The large plusses represent the absolute magnitude limits where 90% 
of companions can be detected at a given physical separation. The numbers indicate the positions of the companions 
listed in Table 5: (1) GJ 84B, (2) GJ 65B, (3) GJ 661B, (4) GJ 257B, (5) GJ 1116B, (6) GJ 860B, (7) GJ 1245B, 
(8) GJ 896B, (9) GJ 1230B, (10) GJ 229B, (11) GJ 618B, and (12) LP 771-95B. The large blank space in the center 
and right-hand-side of the diagram is a clear representation of the "brown dwarf desert", (b) Same as (a), but using 
absolute F110W, and omitting separations < l'.'O. While the sensitivity to T dwarfs is increased in (b), the sensitivity 
to L dwarfs is decreased and close separations cannot be probed. See §5.3 for discussion. 
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Fig. 8. — Statistical semi-major axis distribution for companion search around all 141 M 
dwarf components within 10 pc. The shaded area indicates the separation ranges we consider 
when calculating the multiplicity fraction, with the dashed lines indicating the inner radius 
limits for M and L dwarfs (5 AU), and for the T dwarfs (10 AU) and the outer radius for 
both (70 AU). Search completeness diminishes with increasing separation because NIC2's 
field of view limits our search radius to ~ 9". 
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Table I. Sample Tallies 



Sub-sample Tally 



HST Visits 217 

Total Resolved Targets a 233 

Total Unresolved Targets' 3 22 

Total Star Systems 201 

Resolved Targets within 10 pc a 218 

Unresolved Targets within 10 pc b 21 

Star Systems within 10 pc c 188 

Resolved M Dwarfs within 10 pc a 141 

Systems with M Dwarf Primary within 10 pc c 126 

Resolved Targets beyond 10 pc a 15 

Unresolved Targets beyond 10 pc b 1 

Star Systems beyond 10 pc c 13 



a Each PSF not known to be a blend of more than one 
star or brown dwarf and that is not a background source is 
counted as one "resolved target" . 

b Each PSF known to be the blend of two or more stars 
or brown dwarfs in the same system is counted as one "un- 
resolved target". 

c Denotes known physical association at any physical sep- 
aration, including systems comprising multiple fields of 
view. 
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Fig. 9. — Probability density distributions for select multiplicity fractions listed in Table 6, 
calculated using the binomial distribution. The shaded areas correspond to 68% of the area 
under the curve, equivalent to the la confidence range. The individual plots correspond to: 
(a) M dwarf companions, (b) L dwarf companions, (c) L and T dwarf companions. 
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Fig. 10. — Mass distribution for the binaries in Table 5. The horizontal dashed lines denote 
the hydrogen burning limit (0.075 M ) and the 90% detection limits for this search assuming 
brown dwarf ages of 1 Gyr and 3 Gyr (Table 7). As the masses of the primary components 
approach the hydrogen burning limit, the mass ratios tend to unity, thus implying that 
brown dwarfs rarely form as secondaries. See §6.3—6.4 for discussion. From left to right, the 
binaries are ordered as they appear in Table 5. 
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T4.5 


10.2 


2003 Jan 23 


14.60 


15.6 


17.6 


19.1 


19.6 


20.6 


20.6 


20.6 


20.6 


a, a 


G 99-49 




M3.5V 


5.2 


2003 Jan 8 


6.33 


10.7 


12.2 


12.7 


14.2 


15.2 


16.7 


17.7 


18.7 




LHS 1805 




M3.5V 


7.5 


2003 Mar 8 


6.91 


10.5 


12.5 


13.5 


14.0 


15.0 


16.5 


18.5 


18.5 




LHS 1809 




M5.0V 


9.2 


1998 Feb 12 


8.73 


11.9 


12.4 


14.9 


15.9 


16.4 


17.9 


18.9 


18.9 




GJ 226 




M2.5V 


9.3 


1998 Feb 15 


6.26 


9.4 


10.4 


11.9 


12.9 


13.4 


15.4 


17.4 


17.4 




GJ 229A 


A/B 


M1.5V 


5.7 


1997 Aug 15 


4.19 


8.4 


9.4 


10.4 


11.4 


12.9 


15.4 


17.4 


17.4 


b 


GJ 232 




M4.0V 


8.3 


1997 Dec 20 


8.19 


11.6 


13.6 


14.6 


15.6 


16.1 


18.6 


18.6 


18.6 




L 032-009(A) 


A 


M2.5V 


9.0 


2003 May 7 


5.72 


9.4 


10.4 


11.9 


12.9 


13.4 


15.9 


16.9 


16.9 


b 


L 032-008(B) 


B 


M3.0V 


9.0 


2002 Sep 23 


6.49 


9.7 


10.7 


12.7 


13.7 


14.2 


16.7 


17.7 


17.7 


b 


GJ 244AB 


A/B 


A1.0V 


2.6 


2003 Mar 27 


-1.39 










6.0 


8.5 


9.5 


9.5 


b,d,f,i 


GJ 250B 


B 


M2.0V 


8.7 


1998 Mar 24 


5.89 


8.2 


10.2 


12.2 


13.2 


14.2 


16.2 


17.2 


17.2 


b 
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Name 


Resolved PSFs 


Spectral 


Distance 


Epoch 


App. Mag. 






Absolute 


F180M 


Magnitude Limit 










in FOV 1 


Tvpe 


(pc) 




(F180M) 


0'.'2 


0'.'4 


0'/6 


0'.'8 


l'.'O 


2'/0 


3'.'0 


4'.'0 


Notes 


GJ 257A 


A/B 


M3.0V 


8.0 


1998 Oct 30 


7.15 


10.5 


10.5 


12.5 


13.5 


14.5 


15.5 


17.5 


17.5 


b,h 


GJ 257B 


A/B 


M3.0V 


8.0 


1998 Oct 30 


7.18 


10.5 


11.5 


12.5 


13.5 


14.5 


15.5 


17.5 


17.5 


b,h 


GJ 1093 




M5.0V 


7.7 


1997 Aug 29 


8.56 


12.1 


13.1 


15.1 


16.1 


16.6 


18.1 


18.1 


18.1 




LHS 224AB 


A/B 


M4.5VJ 


9.2 


2003 Mar 13 


8.71 




12.9 


15.4 


15.9 


17.9 


18.9 


19.9 


19.9 


b,c 


GJ 280A 


A 


F5.0IV-V 


3.5 


2003 Jan 11 


-0.67 










6.8 


8.3 


10.3 


11.3 


a,b,d.f ,i 


GJ 283A 


A 


DZQ6 


9.1 


2003 Mar 20 


12.64 


13.8 


15.8 


16.8 


17.8 


18.8 


19.8 


19.8 


19.8 


a,fo 


GJ 283B 


B 


M6.5V 


9.1 


2003 Mar 25 


9.74 


12.9 


14.9 


15.4 


15.9 


16.9 


18.9 


19.9 


19.9 


b 


GJ 1103 




M4.5V 


8.7 


2002 Sep 9 


8.03 


12.3 


13.3 


14.8 


14.8 


15.3 


17.3 


18.3 


18.3 




GJ 293 




DQ9 


7.9 


2002 Aug 29 


12.49 


15.5 


16.0 


17.0 


18.0 


19.5 


21.0 


21.0 


21.0 


a 


GJ 1105 




M4.0V 


8.2 


2003 Mar 10 


7.14 


11.1 


12.1 


13.1 


14.1 


15.1 


16.6 


17.6 


17.6 




GJ 2066 




M4.0V 


6.8 


1998 Sep 17 


5.86 


9.2 


10.2 


11.7 


13.2 


13.7 


14.7 


16.7 


16.7 




GJ 1111 




M6.0V 


3.6 


2003 Mar 2 


7.68 


13.4 


14.9 


15.9 


16.4 


17.4 


18.4 


19.9 


19.9 




GJ 318 




DA6 


8.8 


2003 Jul 4 


11.60 


13.9 


14.9 


15.9 


16.9 


18.4 


21.9 


20.9 


20.9 


a 


GJ 1116A 


A/B 


M5.5VJ 


5.2 


1998 Nov 9 


7.83 


10.2 


12.2 


14.7 


15.7 


16.2 


18.2 


19.2 


19.2 


b,h 


GJ 1116B 




M5.5VJ 


5.2 


1998 Nov 9 


8.17 


10.6 


12.6 


13.1 


14.6 


16.1 


17.6 


18.6 


18.6 


b,h 


LHS 2090 




M6.0V 


6.3 


2003 Jan 25 


8.77 


12.3 


13.8 


15.8 


17.3 


17.8 


19.8 


20.8 


20.8 




GJ 338A 


A 


M0.0V 


6.1 


2002 Nov 27 


4.04 


8.6 


10.6 


11.1 


11.6 


12.6 


14.1 


16.1 


16.1 


b,d,i 


GJ 341 




MO.OV 


10.4 


1997 Oct 15 


5.73 


7.6 


9.6 


11.6 


12.6 


13.1 


14.6 


16.6 


16.6 


g 


GJ 357 




M2.0V 


9.0 


2003 Feb 13 


6.72 


10.9 


12.9 


13.4 


14.4 


14.9 


16.9 


17.9 


17.9 




GJ 1128 




M4.0V 


6.5 


1998 Nov 1 


7.31 


11.7 


12.7 


14.2 


14.7 


16.2 


17.2 


19.2 


19.2 




GJ 367 




M2.0V 


9.7 


1997 Sep 17 


5.99 


8.5 


11.0 


12.0 


13.0 


13.5 


16.0 


17.0 


17.0 




GJ 370 




K6V 


11.1 


1997 Aug 12 


4.86 


8.6 


9.6 


10.6 


11.6 


12.1 


14.6 


15.6 


16.6 


g 


LHS 2206 




M4.0V 


9.2 


2003 Feb 19 


8.66 


10.8 


12.3 


13.8 


15.8 


16.3 


18.8 


19.8 


19.8 




GJ 380 




K7.0V 


4.8 


2003 Feb 19 


3.30 


6.6 


8.1 


9.6 


11.1 


12.1 


14.6 


15.6 


15.6 


d,i 


GJ 388 




M2.5V 


4.8 


1998 Mar 26 


4.81 


8.9 


11.4 


12.4 


13.4 


13.9 


16.4 


17.4 


17.4 




GJ 393 




M2.0V 


7.1 


2003 Jun 16 


5.46 


9.7 


11.2 


12.7 


13.2 


13.7 


16.2 


17.2 


17.2 




LHS 288 




M5.5V 


4.7 


1997 Jul 31 


6.14 


9.2 


11.7 


13.7 


14.2 


15.2 


16.7 


17.7 


17.7 




LHS 292 




M6.5V 


4.5 


2004 Jun 2 


8.26 


11.2 


13.2 


14.7 


15.7 


16.7 


19.2 


20.7 


20.7 


a,d 


GJ 1138AB 


A/B 


M4.5VJ 


9.7 


2003 May 10 


8.04 


10.1 


12.1 


14.1 


15.1 


15.1 


16.1 


17.1 


17.1 


b,c,d 


GJ 402 




M4.0V 


6.8 


1998 Mar 15 


6.67 


10.5 


12.0 


13.5 


14.0 


15.0 


16.5 


17.5 


17.5 




GJ 406 




M5.5V 


2.3 


2003 Feb 27 


6.41 


14.0 


14.5 


16.0 


16.5 


17.5 


19.5 


20.5 


20.5 




GJ 408 




M2.5V 


6.7 


2003 Feb 6 


5.60 


10.0 


11.0 


12.5 


13.5 


14.5 


15.5 


17.5 


17.5 




GJ 411 




M2.0V 


2.5 


2003 Jun 22 


3.72 


10.0 


11.0 


13.0 


13.5 


14.5 


16.0 


17.0 


18.0 


d,i 


GJ 412A 


A 


M1.0V 


4.8 


2003 Jun 14 


5.00 


10.6 


11.1 


12.6 


13.6 


14.6 


15.6 


17.6 


17.6 


b,d,i 



Table 2 — Continued 



Name 


Resolved PSFs 


Spectral 


Distance 


Epoch 


App. Mag. 






Absolute F180M Magnitude Limit 










in FOV 1 


Type 


(pc) 




(F180M) 


0'.'2 


0^4 


0'.'6 


0"8 


l'.'O 


2'.'0 


3('0 


4'.'0 


Notes 


GJ 412B 


B 


M5.5V 


4.8 


2003 May 5 


8.23 


11.8 


14.3 


15.8 


16.8 


17.3 


18.8 


19.8 


19.8 


b 


GJ 432A 


A 


K0.0V 


9.5 


2003 Mar 2 


4.14 




9.1 


10.1 


10.6 


12.1 


13.1 


14.1 


14.1 


h,d,i 


GJ 432B 


B 


DC 


9.5 


2004 Jun 29 


13.67 


15.1 


15.6 


17.1 


17.1 


17.6 


17.6 


17.6 


17.6 


&,b,h 


GJ 433 




M2.0V 


8.9 


2003 Feb 3 


5.76 


8.5 


10.5 


12.0 


13.0 


14.0 


15.0 


16.0 


16.0 




GJ 434 




G8.0V 


9.6 


1997 Dec 13 


3.64 


4.6 


7.6 


10.1 


10.6 


11.6 


13.1 


14.1 


15.1 




GJ 438 




M1.0V 


10.9 


2002 Sep 3 


6.58 


8.8 


10.3 


11.8 


12.3 


13.3 


14.8 


15.8 


15.8 


d,g,i 


GJ 440 




DQ6 


4.6 


1997 Jul 29 


11.20 


16.7 


17.7 


18.7 


18.7 


18.7 


19.7 


19.7 


19.7 


e 


GJ 442A 


A 


G2.0V 


9.2 


1998 Aug 16 


3.30 


5.7 


8.2 


9.2 


10.2 


11.7 


12.7 


14.2 


15.2 


b 


GJ 442B 


B 


M4.0V: 


9.2 


2002 Aug 31 


8.24 


9.7 


11.7 


13.7 


14.2 


15.2 


17.2 


17.2 


17.2 


b 


GJ 445 




M3.5V 


5.3 


2003 Mar 28 


6.24 


9.9 


11.9 


12.9 


13.9 


14.9 


16.4 


18.4 


18.4 




GJ 447 




M4.0V 


3.3 


1997 Jul 13 


5.93 


10.9 


12.9 


14.4 


15.4 


16.9 


18.4 


19.4 


19.4 




GJ 1151 




M4.5V 


8.1 


1997 Jul 15 


7.93 


10.4 


13.4 


14.4 


15.4 


16.4 


17.4 


17.4 


17.4 




GJ 450 




M1.5V 


8.6 


2003 Apr 2 


5.74 


8.8 


9.8 


11.3 


12.3 


13.3 


15.3 


16.3 


16.3 




GJ 451 




K1.0VI 


9.0 


2003 Mar 5 


4.50 


7.7 


9.2 


10.7 


11.2 


12.2 


14.2 


15.2 


15.2 




GJ 1154 




M4.5V 


8.3 


1998 Mar 15 


7.84 


9.9 


12.4 


13.9 


15.4 


16.4 


17.4 


18.4 


18.4 


• ■ -GO 


GJ 475 




GO.OV 


8.4 


1997 Aug 6 


2.80 


4.4 


6.4 


8.9 


9.4 


11.4 


12.4 


14.4 


15.4 


. . 


GJ 479 




M2.5V 


9.6 


2003 Jul 3 


6.30 


9.1 


11.1 


12.1 


12.1 


14.1 


16.1 


17.1 


17.1 


... 


LHS 337 




M4.0V 


6.3 


1998 Sep 3 


7.70 


10.5 


12.5 


14.0 


15.0 


16.5 


17.0 


18.0 


18.0 




GJ 480.1 




M3.0V 


7.9 


2003 Jan 7 


7.68 


9.5 


12.5 


14.0 


15.0 


16.0 


16.5 


16.5 


17.5 




GJ 486 




M4.0V 


8.3 


2003 May 16 


6.66 


9.4 


11.4 


12.4 


12.9 


14.4 


15.4 


17.4 


17.4 




GJ 493.1 




M4.5V 


8.1 


1997 Aug 14 


7.97 


11.0 


13.5 


14.5 


15.0 


16.5 


17.5 


17.5 


17.5 




GJ 494 




MO.OV 


11.4 


1997 Aug 13 


5.74 


8.2 


10.7 


11.7 


12.2 


13.7 


14.7 


15.7 


16.7 


g 


GJ 506 




G7.0V 


8.5 


2003 Jul 1 


2.97 




7.3 


8.3 


9.3 


10.3 


12.3 


13.3 


14.3 


d,i 


GJ 518 




DZ9 


8.2 


1998 Jan 3 


12.81 


14.9 


15.4 


16.9 


17.4 


17.4 


17.9 






a.e 


LHS 2784 




M3.5V 


9.2 


2002 Dec 5 


7.29 


9.2 


11.2 


13.2 


14.2 


15.7 


16.7 


18.2 


18.2 




GJ 551 


C 


M5.0V 


1.3 


2003 Apr 6 


4.84 


10.4 


12.9 


14.4 


15.4 


16.4 


17.4 


19.4 


19.4 


h,d,i 


LHS 2930 




M6.5V 


9.6 


1997 Oct 18 


10.21 


12.1 


15.1 


16.1 


16.6 


17.1 


17.6 


17.6 


17.6 


e 


GJ 555 




M4.0V 


6.2 


1998 Mar 20 


6.19 


10.5 


11.5 


13.0 


13.5 


15.0 


16.0 


18.0 


18.0 




GJ 559A 


A 


G2.0V 


1.3 


1998 Oct 22 


-1.89 








7.9 


8.9 


9.9 


12.4 


12.4 


a.,b,d,i 


GJ 559B 


B 


KOV 


1.3 


1998 Oct 19 










8.4 


9.9 


11.4 


13.4 


13.4 


a,6,z 


G 239-25AB 


A/B 


M3VJ 


9.8 


1998 Nov 7 


6.73 


9.0 


11.0 


12.0 


13.5 


14.0 


16.0 


17.0 


17.0 


b,f,g 


GJ 566A 


A/B 


G7.0V 


6.7 


1998 Sep 14 


2.94 






8.8 


10.8 


12.8 


13.8 


14.8 


14.8 


h,h,i 


GJ 566B 


A/B 


K4V 


6.7 


1998 Sep 14 






6.8 


8.8 


8.8 


10.3 


11.8 






a.,b,e,h,i 


TVLM 513-46546 




M9.0V 


10.5 


2002 Sep 8 


11.16 


11.4 


13.4 


14.9 


15.4 


16.4 


17.4 


17.9 


17.9 
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Name 


Resolved PSFs 


Spectral 


Distance 


Epoch 


App. Mag. 






Absolute 


F180M 


Magnitude Limit 










in FOV 1 


Tvpe 






(F180M) 


0'.'2 


0'.'4 


0'.'6 


0'.'8 


l'.'O 


2'.'0 


3'.'0 


4'.'0 


Notes 


GJ 581 




M3.0V 


6.3 


1998 May 6 


6.07 


9.5 


11.5 


13.0 


14.0 


15.5 


17.0 


18.0 


18.0 




GJ 588 




M2.5V 


5.9 


1997 Sep 18 


4.95 


7.6 


10.6 


12.1 


13.1 


14.6 


16.1 


17.1 


17.1 




GJ 609 




M3.5V 


9.9 


1997 Dec 31 


7.62 


10.5 


13.0 


13.5 


14.5 


15.0 


17.0 


17.5 


17.5 




GJ 618B 


A/B 


M4.5V 


8.3 


1998 Oct 14 


8.83 


9.9 


11.4 


12.4 


12.9 


13.9 


16.4 


16.4 


16.4 


b,e 


GJ 623AB 


AB 


M2.5VJ 


8.0 


1998 Sep 11 


6.16 




10.5 


11.5 


13.0 


13.5 


15.5 


16.5 


16.5 


b,c 


GJ 625 




M1.5V 


6.5 


2002 Sep 10 


5.94 


7.9 


10.9 


11.9 


12.9 


13.9 


15.9 


16.9 


17.9 




GJ 628 




M3.5V 


4.2 


1997 Aug 12 


5.37 


9.8 


10.8 


12.8 


14.3 


14.8 


16.8 


17.8 


17.8 




GJ 631 




K0.0V 


9.7 


1998 Jul 13 


3.69 


7.1 


9.1 


10.1 


11.1 


11.6 


13.1 


14.1 


15.1 




GJ 633 




M2.5V 


21.9 


1997 Sep 6 


8.33 


10.3 


11.3 


12.8 


13.3 


14.8 


15.3 


15.3 


15.3 


g 


GJ 638 




K7.0V 


9.8 


2002 Aug 28 


4.88 


7.0 


9.0 


10.0 


10.5 


11.5 


14.0 


15.0 


16.0 


d 


GJ 643 




M3.0V 


6.4 


2002 Oct 15 


7.03 


11.5 


12.5 


14.0 


15.0 


16.0 


18.0 


19.0 


19.0 




GJ 644ABD 


A/BD 


M2.5VJ 


6.4 


1998 Oct 15 


4.78 




8.0 


10.0 


10.5 


11.5 


13.0 


15.0 


16.0 


h,c,d 


GJ 644C 


C 


M7.0V 


6.4 


2003 Jun 8 


9.20 


12.0 


14.0 


16.0 


16.5 


17.0 


18.0 


18.0 


18.0 


h,d 


GJ 1207 




M3.5V 


8.6 


1997 Sep 16 


7.29 


10.3 


12.3 


13.8 


14.3 


15.8 


17.3 


17.3 


17.3 




GJ 649 




M0.5V 


10.2 


1997 Sep 23 


5.69 


7.9 


9.9 


10.9 


12.9 


14.4 


15.9 


16.9 


16.9 


g 


LHS 3262 




M5.0V 


9.4 


1997 Aug 17 


8.07 


10.1 


13.1 


14.6 


15.1 


16.1 


17.1 


17.6 


17.6 




G 203-47AB 


AB 


M3.5VJ 


7.4 


2002 Nov 8 


6.76 


10.1 


12.1 


13.1 


14.1 


14.6 


16.6 


17.6 


17.6 


b,c,d 


GJ 661AB 


A/B 


M3.0VJ 


6.4 


2002 Aug 26 


5.07 


7.5 


9.5 


11.0 


13.0 


13.5 


16.0 


17.5 


18.0 


b,d,f,h 


GJ 664 




K5.0V 


5.9 


2002 Oct 22 








10.6 


11.1 


12.6 


14.1 


15.1 


15.1 


d 


GJ 666B 


A/B 


K7.0V 


8.7 


1998 Oct 26 


4.95 


8.3 


9.8 


11.3 


12.3 


13.3 


15.3 


16.3 


16.3 


b 


GJ 667A 


A/B 


K4.0VJ 


7.2 


2003 Feb 1 


3.23 


6.7 


7.7 


8.7 


9.7 


10.7 


12.7 


14.7 


15.7 


h,d,h,i 


GJ 667B 


A/B 


K4.0VJ 


7.2 


2003 Feb 1 


3.23 


6.7 


8.7 


9.7 


9.7 


11.7 


13.7 


15.7 


15.7 


h,d,h,i 


GJ 673 




K7.0V 


7.7 


1997 Oct 27 


4.20 


6.6 


8.6 


10.6 


11.6 


13.1 


14.6 


15.6 


16.6 




GJ 674 




M2.5V 


4.5 


1997 Sep 7 


5.02 


9.7 


10.7 


12.7 


13.7 


14.7 


16.7 


17.7 


17.7 




GJ 678.1 




M0.5V 


9.9 


2002 Aug 22 


5.65 


7.0 


10.0 


11.5 


12.0 


13.0 


15.0 


16.0 


16.0 


d 


GJ 682 




M4.0V 


5.0 


1998 Oct 30 


5.79 


8.5 


11.5 


13.5 


14.0 


15.0 


16.5 


17.5 


17.5 




GJ 687 




M3.0V 


4.5 


1998 Feb 5 


4.61 


8.7 


10.2 


12.7 


13.2 


13.7 


15.7 


17.7 


17.7 




GJ 686 




M0.5V 


8.0 


1997 Aug 1 


5.69 


8.5 


10.5 


12.5 


14.0 


14.5 


15.5 


16.0 


16.5 




GJ 694 




M3.0V 


9.5 


1998 Jul 24 


6.10 


9.1 


10.1 


11.6 


13.1 


14.1 


16.1 


17.1 


17.1 




GJ 2130BC 




M2.0V 


14.1 


1998 Sep 17 


6.75 


8.7 


10.2 


11.2 


12.2 


13.2 


16.2 


16.2 


16.2 


b,c,g 


GJ 1221 




DXP9 


6.0 


1997 Aug 15 


12.48 


16.1 


17.1 


18.1 


19.1 


19.1 


19.1 


19.1 


19.1 




GJ 699 




M3.5V 


1.8 


1997 Sep 21 


4.82 


10.7 


12.7 


14.7 


15.7 


17.2 


17.7 


19.7 


19.7 




GJ 701 




M1.0V 


7.7 


1997 Aug 14 


5.45 


8.5 


10.5 


12.5 


13.0 


14.0 


14.5 


15.5 


15.5 




GJ 702A 


A/B 


KO.OV 


5.1 


2002 Sep 18 


1.88 






8.5 


9.5 


10.5 


12.5 


14.5 


14.5 


b,d,i 



Table 2 — Continued 



Name 


Resolved PSFs 


Spectral 


Distance 


Epoch 


App. Mag. 






Absolute F180M Magnitude Limit 










in FOV 1 


Type 


(PC) 




(F180M) 


<y/2 


0'.'4 


0"6 


0'.'8 


l'.'O 


2'/0 


3'.'0 


4'/0 


Notes 


GJ 702B 


A/B 


K5.0V 


5.1 


2002 Sep 18 


1.88 




6.5 


8.5 


10.5 


12.5 


13.5 


14.5 


14.5 


h,d,e,i 


GJ 1224 




M4.0V 


7.5 


1997 Sep 19 


7.98 


11.6 


12.6 


14.1 


14.6 


16.6 


17.6 


18.6 


18.6 




LHS 3376 




M4.5V 


7.2 


2002 Aug 20 


8.26 


12.2 


13.7 


14.2 


14.7 


16.2 


17.2 


17.7 


17.7 




GJ 713AB 


AB 


F7VJ 


8.1 


2003 Feb 28 


2.37 






7.5 


8.5 


9.5 


11.5 


13.5 


13.5 


h,c,d,i 


GJ 1227 




M4.5V 


8.2 


1997 Jul 30 


7.93 


12.4 


13.4 


14.4 


15.4 


15.9 


17.4 


17.4 


17.4 




GJ 721 




A0.0V 


7.7 


2002 Dec 12 


-0.03 








4.1 


5.6 


8.6 


9.6 


9.6 


a,d,i 


GJ 1230AC 


AC/B 


M4.0VJ 


8.2 


2003 Jan 23 


6.91 


8.4 


9.9 


12.4 


12.9 


13.4 


16.4 


17.4 


17.4 


b,c,d 


GJ 1230B 


AC/B 


M5.0V 


8.2 


2003 Jan 23 


8.03 


11.4 


13.4 


14.4 


15.4 


16.4 


17.4 


17.4 


17.4 


h,d,e 


GJ 725A 


A 


M3.0V 


3.5 


1998 Jul 21 


4.59 


9.3 


10.3 


12.3 


13.3 


15.3 


16.3 


17.3 


17.3 


b 


GJ 725B 


B 


M3.5V 


3.5 


1998 Oct 13 


5.15 


8.3 


11.3 


13.3 


14.3 


15.3 


17.3 


18.3 


18.3 


b 


GJ 745A 


A 


K7.0V 


8.6 


1998 Aug 25 


6.68 


9.8 


11.3 


12.3 


13.3 


14.3 


15.3 


17.3 


17.3 


b 


GJ 745B 


B 


M1.0V 


8.6 


2002 Oct 7 


6.75 


9.3 


10.3 


12.3 


13.3 


14.3 


15.3 


17.3 


17.3 


b,d,e 


GJ 747AB 


AB 


M3.5VJ 


8.1 


2002 Dec 11 


6.66 






10.4 


11.4 


12.4 


14.4 


15.4 


17.4 


a,b,c,d 


GJ 752A 


A 


M2.5V 


5.8 


1998 Nov 10 


4.73 


8.7 


10.2 


11.7 


12.2 


13.2 


15.2 


16.2 


17.2 


b 


GJ 752B 


B 


M8.0V 


5.8 


2002 Aug 6 


9.18 


12.2 


14.2 


15.2 


16.2 


17.2 


19.2 


19.2 


19.2 


b 


GJ 1235 




M4.0V 


9.9 


1997 Aug 16 


8.14 


11.0 


13.0 


14.0 


15.5 


16.0 


17.0 


18.0 


18.0 




GJ 764 




G9.0V 


5.7 


2003 Jan 1 


3.04 




8.2 


9.2 


10.2 


11.7 


13.2 


14.2 


15.2 


1 


GJ 768 




A7.0V 


5.1 


1997 Oct 13 


0.18 






8.5 


9.5 


10.5 


12.5 


13.0 


13.5 


1 


GJ 1245 AC 


A/B/C 


M5.5VJ 


4.5 


1998 Oct 10 


7.32 


11.7 


12.7 


14.2 


15.2 


16.2 


17.7 


18.7 


18.7 


b,f,h 


GJ 1245B 


A/B/C 


M6.0V 


4.5 


2003 Feb 14 


7.67 


11.2 


12.7 


13.7 


14.7 


15.7 


17.7 


18.7 


18.7 


b 


GJ 780 




G8.0IV 


6.1 


1997 Aug 25 


1.92 




7.1 


8.6 


10.1 


11.1 


13.1 


14.1 


14.1 


1 


GJ 783A 


A/B 


K2.5V 


6.0 


2003 Mar 21 


3.00 




8.1 


9.1 


9.6 


10.6 


13.1 


14.1 


14.1 


b,d,i 


GJ 783B 


A/B 


M2.5V 


6.0 


2003 Mar 21 




9.1 


11.1 


12.1 


13.1 


14.1 


16.1 


17.1 


17.1 


h e 


GJ 784 




MO.OV 


6.2 


1997 Oct 14 


4.39 


7.0 


10.0 


11.5 


12.0 


13.0 


15.0 


16.0 


16.0 




GJ 785 




K2.0V 


8.9 


1997 Oct 19 


3.73 




8.3 


9.3 


10.3 


11.3 


13.3 


14.3 


14.3 




GJ 1253 




M5.0V 


9.5 


1998 Jan 7 


8.27 


10.1 


12.1 


13.6 


15.1 


16.1 


17.1 


18.1 


18.1 




GJ 791. 2AB 


AB 


M4.5VJ 


8.8 


2002 Oct 28 


7.67 


10.3 


12.3 


13.3 


14.3 


15.3 


17.3 


18.3 


18.3 


b,c,d 


GJ 793 




M3.0V 


8.0 


1997 Aug 9 


6.03 


8.5 


11.0 


12.5 


13.5 


15.5 


16.5 


17.5 


17.5 




GJ 809 




MO.OV 


7.0 


2003 Jul 2 


4.92 


7.8 


9.8 


10.8 


11.8 


12.8 


14.8 


15.8 


16.8 


d,i 


GJ 820A 


A 


K5.0V 


3.5 


2002 Dec 3 


2.54 






9.3 


10.3 


11.3 


13.3 


14.3 


15.3 


h,d,i 


GJ 820B 


B 


K7.0V 


3.5 


2002 Oct 16 


2.89 


7.3 


8.3 


9.3 


11.3 


12.3 


14.3 


15.3 


15.3 


h,d,i 


GJ 827 




F9.0V 


9.2 


1997 Aug 30 


2.96 


6.2 


7.2 


9.2 


10.2 


11.2 


13.2 


14.2 


14.2 


i 


GJ 829AB 


A/B 


M3.0VJ 


6.7 


2002 Nov 17 


5.74 


9.4 


10.9 


11.9 


12.9 


13.9 


15.9 


16.9 


16.9 


b,c,d 


GJ 831AB 


A/B 


M4.0VJ 


7.9 


2002 Oct 11 


6.69 


8.5 


10.5 


12.5 


13.5 


14.5 


16.5 


17.5 


17.5 


b,c 



Table 2 — Continued 



Name 


Resolved PSFs 


Spectral 


Distance 


Epoch 


App. Mag. 






Absolute F180M Magnitude Limit 










in FOV 1 


Type 


(pc) 




(F180M) 


0'.'2 


0"4 


n"fi 


U. o 


i"n 
1 . u 


9"n 


3'.'0 


4'.'0 


Notes 


GJ 832 




M1.5V 


4.9 


1997 Jul 28 


4.60 


8.5 


10.5 


11.5 


12.5 


14.5 


15.5 


16.5 


16.5 




G 188-38 




M3.5V 


8.9 


1997 Oct 11 


6.96 


8.7 


10.2 


12.2 


13.2 


14.2 


15.2 


16.2 


16.2 




GJ 846 




M3.5V 


10.2 


1997 Nov 17 


5.36 


8.0 


9.0 


11.0 


12.0 


13.0 


14.0 


15.0 


16.0 


§ 


LHS 3746 




M3.5V 


7.4 


2002 Nov 18 


6.79 


9.6 


10.6 


11.6 


13.6 


14.6 


16.6 


17.6 


17.6 




GJ 845A 


A 


K4.0V 


3.6 


1997 Aug 4 


2.17 






10.2 


11.2 


12.2 


14.2 


16.2 


16.2 




GJ 849 




M3.0V 


8.9 


1997 Nov 15 


5.69 


9.2 


10.2 


11.2 


12.2 


13.2 


15.2 


17.2 


17.2 




LHS 3799 




M4.5V 


7.4 


1997 Oct 29 


7.47 


10.1 


12.6 


14.1 


14.6 


15.6 


16.6 


17.6 


17.6 




GJ 860A 


A/B 


M3.0V 


4.0 


1998 Nov 10 


5.04 


10.0 


11.0 


13.0 


14.0 


15.0 


17.0 


18.0 


18.0 


b d 


GJ 860B 


A/B 


TV T A 7 

M4.0V 


4.0 


1998 Nov 10 


5.04 


12.0 


13.0 


14.0 


14.5 


15.0 


17.0 


18.0 


18.0 


h.d e 


GJ 867AC 


AC 


M2.0VJ 


8.6 


2002 Nov 18 


5.11 


8.8 


10.3 


11.3 


12.3 


13.3 


15.3 


16.3 


16.3 


b c d i 


GJ 867B 


B 


M3.5V 


8.6 


1998 Aug 14 


6.66 


8.8 


10.3 


12.3 


12.8 


14.3 


16.3 


17.3 


17.3 


b 


GJ 873 




M3.5V 


5.0 


2002 Dec 24 


5.55 


8.5 


10.5 


12.5 


13.5 


14.5 


16.5 


17.5 


18.5 




GJ 876 




M3.5V 


4.6 


1997 Dec 5 


5.16 


8.7 


10.7 


12.7 


13.7 


14.7 


15.7 


17.7 


18.7 




GJ 1276 




DZ9+ 


8.5 


1997 Sep 21 


13.56 


15.3 


16.3 


17.3 


17.3 


17.3 


18.3 


18.3 


18.3 




GJ 877 




M2.5V 


8.6 


1997 Aug 29 


5.94 


8.3 


10.3 


11.3 


12.3 


13.3 


16.3 


17.3 


17.3 




GJ 880 




M1.5V 


6.8 


2002 Oct 15 


4.80 


8.3 


9.8 


10.8 


11.8 


12.8 


14.8 


15.8 


16.8 


d i 


GJ 881 


A 


A4.0V 


7.6 


1998 Aug 6 


0.94 






5.6 


6.6 


7.6 


8.6 


10.6 


11.6 


a,,b,d,i 


GJ 884 




K7.0V 


8.2 


1997 Dec 28 


4.44 


7.4 


9.4 


10.4 


11.4 


12.4 


13.4 


15.4 


15.4 




GJ 887 




M1.0V 


3.2 


1997 Sep 11 


3.45 


7.4 


9.9 


11.4 


12.9 


14.4 


15.4 


17.4 


17.4 




GJ 892 




K3.0V 


6.5 


1997 Sep 19 


3.17 




8.4 


9.9 


10.9 


11.9 


12.9 


14.9 


14.9 




GJ 896A 


A/B 


M3.5V 


6.2 


2004 Jun 30 


5.57 


7.0 


8.5 


10.0 


11.0 


12.0 


14.0 


15.0 


16.0 


a,b 


GJ 896B 


A/B 


M4.5V 


6.2 


2004 Jun 30 


6.55 


8.0 


10.0 


12.0 


12.0 


13.0 


14.0 


15.0 


16.0 


a, 6 


GJ 1286 




M5.0V 


7.2 


1997 Sep 18 


8.42 


13.7 


14.7 


15.7 


16.2 


16.7 


18.7 


18.7 


18.7 




GJ 902 




K3V 


11.4 


1997 Oct 19 


4.56 


7.7 


8.7 


10.7 


11.7 


11.7 


13.7 


14.7 


15.7 


g 


GJ 905 




M5.5V 


3.1 


2003 Jan 12 


6.25 


11.5 


13.5 


14.0 


14.5 


15.5 


17.5 


18.5 


19.5 




GJ 1289 




M3.5V 


8.1 


1997 Sep 10 


7.39 


10.5 


12.5 


13.5 


14.5 


15.5 


17.5 


18.5 


18.5 




GJ 908 




M1.0V 


5.9 


1997 Sep 9 


5.13 


9.6 


10.6 


12.1 


13.1 


14.1 


16.1 


17.1 


17.1 





X A/B denotes that components A and B are resolved. AB denotes an unresolved known multiple. 
a No PSF subtraction 

b Known multiple system, excluding planets 

c Sensitivity measured around unresolved or very close binary 



d 2MASS H magnitude for primary 
° Off-center by more than 5" 

'Sensitivity measured around brighter component only 
g Beyond 10 pc 

h Sensitivity measurements exclude 30° between components. 
'Saturated core 
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Table 3. Background Sources with Companion-Like Colors 



Primary 


Separation 


P. A. (deg) 


Mimics 


Deciding Factor 


GJ 633 


8"5 


157.3 


Late L in F207M-F222M 


Too red to be late L in F110W-F180M 


GJ 633 


3"3 


268.6 


Late L in F207M-F222M 


Too red to be late L in F110W-F180M 


GJ 1093 


8"0 


250.6 


Mid M in F110W colors 


Background Ml or earlier in F207M-F222M 


GJ 1224 


12'.'6 


91.4 


Mid L in F110W-F222M 


Background early M in F180M-F222M 


GJ 367 


8"4 


3.3 


Early L in F180M-F207M 


Background main sequence in F110W-F222M 


GJ 438 


13'.'2 


254.0 


Hot white dwarf 


Too red to be hot WD in F110W-F180M 



Table 4. Astrometry of Unsaturated Resolved Systems 





Pair 


P 


<J P 


PA. 


o-p.A. 


Epoch 




A Magnitude 








(") 




(° E of N) 






F110W 1 


F180M 


F207M 


F222M 


GJ 


1005AB 


0.329 


0.008 


234.4 


1.0 


2002.7532 




1.27 


1.31 


1.32 


GJ 


65AB 


1.653 


0.008 


103.3 


0.2 


2002.8540 




0.18 


0.15 


0.16 


GJ 


84AB b 


0.443 


0.006 


103.4 


1.0 


2002.7506 


4.59 


4.01 


4.18 


3.82 


GJ 


105AC 


3.220 


0.036 


293.5 


0.4 


1998.0225 


5.96 


5.87 


5.45 


5.40 


LP 


771-95AB 


7.706 


0.008 


315.0 


0.0 


2003.4620 




0.36 


0.36 


0.35 


LP 


771-95BC 


1.344 


0.008 


138.1 


0.2 


2003.4620 




0.56 


0.54 


0.50 


GJ 


169. 1AB 


9.201 


0.008 


63.7 


0.0 


2002.7693 




5.89 


5.98 


6.07 


GJ 


229AB 


7.627 


0.031 


164.1 


0.1 


1997.6202 


9.63 


11.44 


9.44 


10.75 


GJ 


257AB 


0.560 


0.008 


280.4 


0.6 


1998.8271 




0.03 


0.00 


0.00 


GJ 


1116AB 


1.498 


0.007 


102.7 


0.2 


1998.8562 


0.42 


0.35 


0.32 


0.30 


GJ 


618AB 


5.574 


0.008 


226.8 


0.0 


1998.7837 




2.66 


2.58 


2.49 


GJ 


644A-BD 


0.258 


0.048 


151.3 


7.6 


1998.7868 




0.41 


0.41 


0.43 


GJ 


661AB 


0.647 


0.008 


187.6 


0.5 


2002.6489 




0.43 


0.29 


0.28 


GJ 


1230AC-B 


5.117 


0.008 


6.1 


0.0 


2003.0605 




1.91 


1.88 


1.85 


GJ 


747AB 


0.234 


0.048 


87.4 


8.4 


2002.9441 




0.02 


0.05 


0.03 


GJ 


1245AB 


6.964 


0.007 


82.6 


0.0 


1998.7734 


0.17 


0.16 


0.15 


0.13 


GJ 


1245AC 


0.594 


0.007 


269.6 


0.5 


1998.7734 


1.49 


1.29 


1.18 


1.08 


GJ 


860AB 


3.184 


0.008 


99.5 


0.1 


1998.8590 




1.01 


0.97 


0.94 


GJ 


896AB 


5.351 


0.008 


87.9 


0.0 


2004.4956 




1.54 


1.48 


1.43 



a Missing values correspond to stars with central core saturation for which a PSF fit is not available. 
b Values from G04, Table 4. 
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Table 5. M Dwarf Multiple Systems with Separations Between 5 and 70 AU Recovered in 

the Survey 



A Component a Spectral 




Mass c 


B Component 


Spectral 




Mass c 


Angular 


Inferred 




Type 




M 






Type 




M Q 


Separation 


Semi-major Axis (AU) f 


GJ 229A 


M1.5V 


5.41 


0.59 


GJ 


229B 


T6.0 


16.81 


0.05: d 


7'.'62 


55.2 


GJ 84A 


M2.5V 


5.95 


0.52 


GJ 


84B 


M7.0V 


10.54 


0.084 e 


0'.'44 


5.2 


GJ 661A 


M3.0V 


6.03 


0.51 


GJ 


66 IB 


M3.0V 


6.03 


0.51 


0'.'70 


5.6 


GJ 896A 


M3.5V 


6.60 


0.39 


GJ 


896B 


M4.5V 


7.62 


0.22 


5'.'35 


42.1 


GJ 618A 


M2.5V 


6.62 


0.38 


GJ 


618B 


M4.5V 


8.43 


0.15 


5'.' 62 


58.7 


GJ 860A 


M3.0V 


7.02 


0.31 


GJ 


860B 


M4.0V 


7.95 


0.18 


3'.' 19 


16.2 


GJ 1230A 


M4.0V 


7.34 


0.26 


GJ 


1230B 


M5.0V 


8.46 


0.15 


5" 11 


53.2 


LP 771-95A 


M2.5V 


7.56 


0.23 


LI' 


771-95B 


M3.5V 


7.92 


0.18 


7'.' 74 


67.8 


GJ 257A 


M3.0V 


7.63 


0.22 


GJ 


257B 


M3.0V 


7.66 


0.21 


0'.'57 


5.7 


GJ 1245A 


M5.5V 


9.05 


0.12 


GJ 


1245B 


M6.0V 


9.40 


0.10 


7'.'01 


40.1 


GJ 1116A 


M5.5V 


9.24 


0.11 


GJ 


1116B 


M6.0V 


9.58 


0.10 


1'.'51 


9.9 


GJ 65A 


M5.5V 


9.32 


0.11 


GJ 


65B 


M6.0V 


9.50 


0.10 


1?66 


5.6 



a Ordered by decreasing mass, as shown in Figure 9. 

h H band photometry from 2MASS. Close binaries were deconvolved adopting AH = AF180M . 
c Based on the Mass-Luminosity Relation of iHenrv fc McCarthy! jl993l) . 

Estimate based on lAIIard et all Jl996h . 
e G04 

statistically corrected for projection effects iFischer fc Marcvlll992h . 

Table 6. Summary of M Dwarf Multiplicity Fractions 



Companion Range 
Spectral Type 


Search Radius 
AU 


Systems 
Probed 


Detections 


Hidden 


Mult. Fraction 

% 


Volume 
Limited? 


MOV - M9V 
L0 - L9 
L0 - L9 
L0 - T9 
L0 - T9 


5 - 70 
5 - 70 
12 - 70 
10 - 70 
14 - 70 


126 
51 
126 
43 
126 


11 



1 
1 


2 






W.3± 3 2 i 

o.atli 

0+ 14 

u - u -o.o 

9 q+5.0 
z -°-0.7 

o.8lJ;f 


Yes 
No 
Yes 
No 
Yes 



Table 7. Brown Dwarf Masses (M ) Based on Models of IChabrier et all (120001 1 



T *f a 


L3 
2000 


L5 
1750 


L8 
1500 


T5 
1200 


T7 
900: 


1 Gyr 
3 Gyr 
5 Gyr 


0.070 
0.073 
0.075 


0.060 
0.070 
0.072 


0.050 
0.057 
0.065 


0.050 
0.052 
0.065 


0.030: 
0.040: 
0.050: 



HGolimowski et alj j2004al) 



